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Abstract 


The primary purpose of this publication is to develop a mathematical model 
to describe smooth paths along any combination of circles and tangent lines. Two 
consecutive circles in a path are either tangent ( externally or internally ) or they 
appear on the same (lateral) or opposite (transverse) sides of a connecting tan- 
gent line. A path may start or end on either a segment or circle. The approach is 
to use mathematics common to robotics to design the path as a multilink manipu- 
lator. This approach allows a hierarchical view of the problem and keeps the 
notation manageable. A user simply specifies a few parameters to configure a 
path. Necessary and sufficient conditions automatically ensure the consistency of 
the inputs for a smooth path. Two example runway exit paths are given, and an 
angle to go assists in knowing when to switch from one path element to the next. 


1. Introduction 

A NASA goal is to improve the aviation system 
throughput in all weather conditions. One aspect of 
this problem is guidance along the runway and exit 
path. In this respect, coordinates along the centerline 
of a runway and the exit path must be known. 
Although others (refs. 1, 2, and 3) have modeled spe- 
cific exit paths (taxiways) with linear segments and/or 
circles, a general formulation does not appear to exist 
for such paths. 

This publication completely characterizes smooth 
paths along an arbitrary combination of circles and 
tangent lines. The approach herein is to configure the 
path as a multilink manipulator and use the mathemat- 
ics common to robotics. This approach provides a 
ready-made high-level perspective of a complex sys- 
tem. A user simply specifies a few parameters in a 
computer program to configure a new path. Two 
examples of designing runway exit paths are given. 
Also, as a convenience, this publication introduces a 
parameter to tell when to switch from one path ele- 
ment to the next. 

2. Symbols 

A' j homogeneous transformation matrix 

from axis system / to axis system j or 
transform from circle C ■ to circle C . 

a j component of vector from C -_j to C ( 

along X- 

C circle 


ccw counterclockwise 

cir circle 

cw clockwise 

Cj circle classification, either exterior 

(/ = 1) or inward (/ = -1) 

(hj,kj, 0) center of circle C • in general refer- 

ence frame (X G ,Y G ,Z G ) 

i, j index 

Lg 1 transformation from circle C- +1 to 

general reference frame with 

'/ + 1 = o 

f length of segment P ■ , positive or 

negative 

N number of circles in path 

P ■ tangent segment before circle C • 

p’j 1 vector (in axis system i - 1 ) from axis 

system i — 1 to axis system i or vector 
from center of circle C j to center of 
circle C • 

Rj radius of circle C ; 

r'j Euler rotation matrix about Z ■ from 

axis system i to axis system j 

seg segment 



(X g ,Y g ,Z g ) 

general reference frame 

(*„7„Z,) 

axis system rigidly attached to center 
of circle C • 

(X'jfY'jfZ'j) 

axis system parallel to (X G ,Y G ,Z G ) 
at center of circle C ■ 

x G*y G ,Z G 

coordinates in general reference 
frame ( X G ,Y G ,Z Q ) 

X G, exit ’-^G, exit 

exit coordinates in general reference 
frame ( X G ,Y G ,Z G ) 


coordinates in circle axis system 
(Xi,Y,Zi) 

0) 

coordinates in circle axis system 
(X',7',Z') parallel to (X G ,Y G ,Z G ) 

(x*,y*,o) 

starting point of motion on circle em- 

0) 

ending point of motion on circle C ( 

r 

azimuth of initial motion on first cir- 
cle C () relative to X G or angle 
between X G and Y 0 

V 

rotation angle about center of circle 
C ■ or rotation about Z -_ j to rotate cir- 
cle C ■ about circle C / _ 1 

A ,- 

central angle of circle C - , measured 
between current motion position and 
starting point for motion on circle; 
negative value means prior to starting 
point 

K 

travel angle (central angle) on circle 
C-_| about Zj j 

V 

central angle on C • between X- and 
radial ray to known point or radial 
projection angle 

V 

central angle on C- between T' and 
radial ray to known point or radial 
projection angle 


Vg,i 

azimuth on circle C • relative to axis 
*G 

Vi 

azimuth on circle C • from north 

Vp. 

azimuth on segment P ■ from north 

Vx G 

azimuth of axis X G 

<°i 

index of rotation, clockwise (i = -1) 
or counterclockwise (j = 1) on circle 

c,-- 1 

Subscripts: 


max 

maximum 

target 

target value on path 


3. Circle Classification 

A known path can be approximated by measuring 
a series of points along the path and then linearly con- 
necting the points (line segments). However, this 
requires many points to approximate closely a curvi- 
linear path, and even then, the slope of the path where 
the segments connect is discontinuous. 

The method discussed in this publication con- 
structs a smooth path by using a combination of cir- 
cles and lines. Any two successive circles in a path are 
either tangent or connected by a common tangent line. 
A path may begin or end on a circle or line. A path is 
the composite of pieces of travel along circles and tan- 
gent lines. A user specifies information about the cir- 
cles and segments to generate the path. If two 
immediately consecutive circles in a path are tangent, 
except for the tangent point, either (1) the second cir- 
cle is external to the first (disjoint interiors), (2) the 
second circle is inside the first circle, or (3) the first 
circle is inside the second circle. 

Any two successive circles in a path (say, first and 
second circles) are either tangent or connected by a 
single co mm on tangent line. The second circle is 
either an exterior circle or inward circle with respect 
to the first circle. 
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Figure 1. Possible locations and classifications of second 
circle. C._j is first of two circles; C- is second circle 
choice. 

3.1. Exterior Circle 

If two circles lie on opposite sides of a given tan- 
gent line, the second circle is an external (or 
transverse) circle with respect to the first, as shown in 
figure 1. 

3.2. Inward Circle 

If the first and second circles lie on the same side of 
a given common tangent line, the second circle is an 
inward (or lateral) circle with respect to the first circle, 
as shown in figure 1 . Notice that the second circle can be 
inside the first or vice versa. 

3.3. Parameter for Circle Type 

Symbolically, the circle classification is 

f 1 (Cj exterior) 

' {-1 (Cj inward) 




Figure 2. Multilink manipulator. 

Designing a path is somewhat analogous to con- 
figuring a planar multilink manipulator if one thinks 
of the manipulator body as the path. Imagine, for 
example, that the manipulator in figure 2 is lying flat 
on a horizontal plane and think of the body of the 
manipulator as a road or path on this plane. A path of 
linear segments then follows. However, for a smooth 
path the segments should make a smooth transition 
into one another. Hence, imagine the joints as centers 
of circles and the links as tangent to these circles. 
Now, imagine moving along a path from the base of 
the manipulator toward the shoulder along link 0, 
which ends as a tangent to the shoulder circle. Proceed 
around the circumference of the shoulder circle for so 
many degrees (travel angle) and then travel toward the 
elbow along link 1, which is tangent to both the shoul- 
der circle and the elbow circle, and so forth until the 
end of the path. Changing the joint angles changes the 
path. Really there are no joints on the path, but this 
publication relates the travel angle along the circular 
arcs to the joint angles to use mathematics common to 
robotics, which is the basic idea in this publication. 


4. Modeling Path as Multilink 
Manipulator 

A robot manipulator (ref. 4) is a series of joints 
connected by links; for example, the sketch in figure 2 
shows a manipulator ann with three rotational joints 
(shoulder, elbow, and wrist). Each joint only moves 
the part of the manipulator beyond that joint; that is, 
the elbow joint rotates the forearm but the wrist joint 
does not. To configure a manipulator, one changes the 
joint angles. 


A known path can be approximated with a path 
along circles and tangent lines. To describe a path 
using circles and lines, one must know the radii of the 
circles, where the circles are located, the angular travel 
along each circle, and the length of any linear segment 
between circles. 

4.1. Axis Systems 

Figure 3 shows the initial layout of circles and linear 
segments to form a path. All axis systems are 
right-handed. As shown, all circles are exterior and all 
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Parallel to X G 



*N-l 


Figure 3. Zero-reference configurations. All circle angles are 0. 


segments are positive, but this initial configuration 
changes, depending on the type of circle and whether a 
segment is positive or negative. To change a circle from 
exterior to inward, simply move the external circle to the 
opposite side of its upper tangent line; to change a posi- 
tive segment to a negative segment, simply shift the seg- 
ment by its length in the opposite direction. As shown, 
there is a general reference frame (Xq plus 

each circle has its own rigidly attached axis system (at 
the center of the circle). When T = 0, the axes X G and 
Yq are aligned. Note that the circle indices begin with 
zero; that is, the axis system of the first circle is 
(Xq,Y () ,Zq). All the axis systems for the circles are par- 
allel when the circles angles are 0. The configuration of 


the circles changes to accommodate different path 
designs. 

4.2. Circle Rotation Angle 

Think of the center of each circle as a rotational 
joint and the distance between consecutive circle 
centers as a link (like in a manipulator robot). The 
first circle is stationary. A rotation at the center of the 
first circle (6 l about Z 0 ) translates and rotates the sec- 
ond circle and all remaining circles. If two circles hap- 
pen to be tangent, the second circle moves along the 
circumference of the first circle. By convention, the 
circle angle 0 ^ is always about Z | . See appendix A 
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for an example of configuring a path with circle 
angles. Later, these circle angles will configure a path 
but will not be user inputs. (With respect to manipula- 
tors, circle angles are analogous to the joint angles.) 

4.3. Segment Notation 


Symbolically, denote the direction of motion on a cir- 
cle as 


°>i = 



(C/ , ccw)j 
(C i _y CW)J 


( 2 ) 


If there is a segment prior to the first circle C 0 , call 
it P 0 . Afterwards, segment P. +1 follows circle C- 
Therefore, the last segment P N follows the last circle 
C N _i . Hence, there are N circles in a path and up to 
N + 1 linear segments. 

5. Smooth Path Design 

Figure 4 shows the six possible locations of the 
second of two immediately consecutive circles in a 
path. Various paths can be traced (arrows) along the 
configuration; however, not all are smooth. (A cusp 
can occur at the transition point between two circles or 
between a circle and a segment; otherwise, an abrupt 
change in direction occurs.) Smooth paths (no cusps) 
are the primary interest in this publication, but the 
basic analysis is not limited in this respect. 

5.1. Direction of Motion on Circle 


5.2. Necessary and Sufficient Conditions for 
Smooth Path 

A computer program can automatically assign 
information or issue an alert for invalid inputs in the 
design of a path. For a smooth transition from a circle 
to a segment, the segment length parameter with alge- 
braic sign is 

/ 0 = Absolute length of segment P {) j 
1/ = co ; ( Absolute length of segment P.) J 

O' = 1,2, ...,N) (3) 

For a smooth path, another relationship that must 
hold between immediately successive circles in a path 
is (fig. 4) 


A counterclockwise (ccw) motion on a circle is 
positive, and a clockwise (cw) motion is negative. 


where (Oj detennines the direction of motion on the 
first circle and depends on the initial path direction 


Counterclockwise Clockwise 



Inward 

circles 

t 

T 

Exterior 

circles 


Figure 4. Possible paths (arrows) along next circle. 
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(right or left). Equations (3) and (4) constitute neces- 
sary and sufficient conditions for a smooth path. 


( 6 ) 



As an example in using equation (4), if the motion 
on the first of two consecutive circles is clockwise 
(co. = -1) and a counterclockwise motion (®. + , = 1) 
is desired on the second circle, then the second circle 
must be external, as shown by 



For the same first circle ( co l ■ = - 1 ), if the second circle 
were external (c- +1 = 1), then the direction of motion 
on the second circle must be counterclockwise, as 
shown by 

®/ +1 = ~ e M °>i = -(D(-D = 1 

Consequently, beyond the first circle, a user only 
needs to specify either the circle classification or the 
direction of motion on the remaining circles. 

6. Homogeneous Transformation 
Matrices Along Path 

A homogeneous transformation matrix (ref. 4) 
completely characterizes the geometric relationship 
between two axis systems, which may be rotated 
and/or displaced from each other. 

6.1. Transformation Between Consecutive 
Circles 

The homogeneous transformation matrix from cir- 
cle Cj to the previous circle C-_j (for i = 1, 2, ..., AO is 


is the component of the vector from C -_, to C- along 
Xj. For the special case of tangent circles (/• = 0) , the 
magnitude of a • is the distance between circle centers. 
The transformation matrix between any two arbitrary 
circles in a path is a concatenation of the individual 
transformation matrices between pairs of circles in a 
path. See appendix B for a brief discussion of homo- 
geneous matrices. 

6.2. Transformation From Circle to General 
Reference Frame 

The homogeneous transformation matrix from any 
circle C ; in the path to the reference frame is 

4 - 444 - 4 1 < 7 > 

where, for the particular orientation of the general ref- 
erence frame shown in figure 3, 

cos T -sin T 0 h 0 

A° g = sin r cos r 0 k 0 ( 8 ) 

0 0 10 

0 0 0 1 _ 

is the transformation from the first circle C () and 

h o = x G exit - (Exit direction )R 0 sin T (9) 

*0 = 34 exit + (Exit direction)/^ cos T (10) 


cos dj -sin 0- 0 a i cos 0. -/■ sin 0 
sin 9j cos 9 1 0 a • sin 0. + / ■ cos 0 

0 0 1 0 

0 0 0 1 

where 0. is the circle angle (positive about Z (1 ) 
from 4 ■_[ to X p /• is the distance (either plus or 
minus) along the tangent segment between the two cir- 
cles, and 




are the center coordinates of C 0 in the reference frame. 
Here, x G it and y G it are the coordinates of the 
exit point (beginning point on first circle in the path). 
The exit point is where the path first veers away from 
a reference direction (along runway), and the parame- 
ter that categorizes the exit direction is 



(Exiting to right) 
(Exiting to left) 
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Exit direction = 


( 11 ) 



6.3. Recursive Relationship for Circle Center 

The center of a circle is important in creating a 
path and later in following a path. Let (/?■,£■, 0) be the 
center point of circle C -. An equation for the circle 
center (zero coordinates in circle axis system) is 


where on the first circle C () , R (] is the radius, 
( x G exit , y G exj| ,0) is the initial point (exit), and (» l is 
the direction of travel. 

7. Drawing Path 


h i 


' 

0 

k i 

, i—1 , i 

• = A G A i-l 

> o 

0 


U 

. 1 . 


1 


Several input formats could be used to define a 
path (user-preference issue). In any event, directly or 
indirectly from a user’s input, the following informa- 
tion about a path is known: 

1 . Exit point 


= A 


i—1 


[ a i cos 0- - /■ sin 0 ( 
a ■ sin 0. + /■ cos 0 ( 

0 

1 


( 12 ) 


2. Whether the path is to the right or left 

3. The number of circles 

4. Whether each circle is exterior or inward 


Or, more explicitly, the following vector recursive 
relationship is given: 


h i 

„/-l 

' ' 

a- cos 9 i - lj sin 0 ( 

k i 

• - r g ■ 

a. sin 0. + /■ cos 9- r 

. 0 


0 


h i - 1 

_i— 1 

' ' 

Oj cos 0j - lj sin 0j 

\ k 
K i - 1 

> + R G < 

dj sin 0j + lj cos 0j 

. 0 . 


0 


To start the recursion, 




0 

h o 


A"(0,3) 

< k 

> = < 

o t 

K o 


Ag(l,3) 

. 0 


0 

- j 


(13) 


(14) 


5. Radius of each circle 

6. Distance along each linear segment (with 
direction) 

7. Angular travel (travel angle with direction) 
along each circle 

A path can be drawn once a set of design parame- 
ters is known. Several reasons for drawing a path are 

1. To show users what their design looks like, 
thereby adding confidence (validity) to the 
design 

2. To provide visual feedback in the design pro- 
cess to changes in parameters 

3. To overlay and adjust the path to known geom- 
etry (terrain) 

4. To provide a graphical reference for motion 
(visual deviations from a path) 


or equivalently, 


h o 


X G,ex it 


Rq sin T 

k 

K o 

> = < 

y G,exit 

- fO| 

Rq cos r 

. 0 


0 


0 


7.1. Path Along Circle 

Figure 3 shows the configuration of the circles and 
segments in a path (for exterior circles and positive 
segments) when the circle angles are zero. Changing 
these circle angles from their zero values moves the 
circles (and segments) to other locations to 
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accommodate a particular path design. However, a 
more natural way to design a path is to specify the 
angular travel (travel angle) along each circle in the 
path instead of specifying the circle angles 0-. (Anal- 
ogously, in robotics, it is more natural for an operator 
to command the movements of an end-effector or 
gripper directly than to specify the joint angles indi- 
vidually in a manipulator.) 

7.1.1. Relating Circle Angle to Travel Angle 

Define the travel angle A- as the angular travel 
along the circumference of circle C •_ j (that is, arc 
travel along the circumference of a circle subtends the 
travel angle). A counterclockwise travel angle is 
positive. 


Notice in figure 5 that motion begins from the rear 
(point on circle along -X-_j ) of an exterior circle and 
from the front (point on circle along X • , ) of an 
inward circle. The circle angle 0-, which moves circle 
Cj with respect to the circle C- , , is a function of the 
travel angle and the circle classification. From 
figure 5, this relationship is seen to be 

e, = 90°( I + c ( _! ) + Xj (16) 

Hence, instead of specifying the circle angles to con- 
figure a path, a user specifies the angular travel along 
each circle and then calculates the circle angles, which 
is more convenient. 

7.1.2. Path Coordinates Along Circle 




Figure 5. Geometry relating arc travel and circle angle for 
exterior and inward circles. 


Let the circle C • be either an exterior or inward cir- 
cle that rotates about circle C ■_ j by the angle . The 
path along circle C j _ l is the same whether there is a tan- 
gent line from C-_j to C-. Therefore, for this 
calculation, assume that the two circles are tangent so 
that circle Cj touches and rotates along the circumfer- 
ence of circle C •_ j. As A • varies, a plot of the touch 
points is the path along circle C,_| (fig. 6). In the 
reference frame, the path coordinates along circle C •_ j 
are given by 


x G,i - 1 

. - V G,/— 1 

0 

1 





*r 

i 

.7-1 

' = a g 

V 1 

* 

o 



0 




1 


(Path along C ■_, ) 


(17) 


where the vector on the right-hand side of 
equation (17) is the touch point in the axis system of 

circle C, and the term in brackets varies with .0 

i 

which varies with A •. The segment length /■ in A ._j 
is set equal to 0 in this calculation, which changes at 
most two elements. Thus, equation (17) gives points 
along the circular track as an implicit function of the 
travel angle. To draw the circular track on each circle 
C ■_ j (for i = 1,2, ..., AO, vary the travel angle A- 
from 0 to its final value. If desired, equation (17) 
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X,. 



m 

Always touch 
point witli 
respect to 
tangent circle, 
exterior or 
inward 


Figure 6. Coordinates on circular path which depend on 
travel angle. 


expands to the more explicit vector expression as 
follows: 


In essence, equation (17) uses the next circle to 
draw a path on the current circle. Consequently, a fic- 
titious circle is required to draw the path along the last 
circle C ( i . For i = N, equation (17) should give the 
path along the last circle , , but up to now, there is 
no matrix A N _ j because there is no circle C N . To use 
equation (17) to draw the path along the last circle 
Cyy_j, think of a fictitious circle C N with arbitrary 
radius R N at the end of the final segment P,,. The 
homogeneous transformation for this fictitious circle 
is 

cos dy ~sin 0 V 0 a N cos 9 N - l N sin d N 
sin 9 n cos 9 n 0 a N sin 9 N + /,. cos 9 N 
0 0 1 0 

0 0 0 1 

(19) 

where 



x G,i - 1 


J5- 

7' 

y G,i-l 

• = < 

k 

K i~ i 

. 0 


. 0 . 


x 


-e -Rf cos 0. + a- cos 6- 
-e/Rj sin 0 / + sin 0 ( 


a N ~ R N-\ + e N R N 
~ 90°(1 + e N _ x ) + 

Here, A N is the angular travel on circle C N _ j, and 
0 N is the angular rotation angle of circle C N about 
circle C N _ j. Computation of this transformation is 
transparent and automatic with respect to a user. 



h i - 1 


Rj_ j cos 0/ 


1 


\ k 
K i - 1 

+ r g ' 

R/_l sin 0 j 

, 0 


0 


which means that a vector from the reference frame to 
a path point on circle C-_ 1 is equal to the vector from 
the reference frame to the center of circle C -_j plus 
the vector from the center of circle C -_ j to the path 
point (where the matrix multiplication transfonns the 
components of the latter vector to the reference 
frame). The initial point on a circular track is where 
the travel angle is 0, and the final point is where the 
travel angle takes on its final value. If the segment 
length between two circles is 0 then the ending point 
on the first circle is the starting point on the second 
circle. (See appendix C for specific equations to 
compute the beginning or ending points on any circle 
in a path.) 


7.2. Path Along Segment 

A straight-line segment can occur 

1 . Prior to the first circle 

2. Between two circles 

3. After the last circle 

In the first case, the initial point is known. Connecting 
this initial point and the exit point (which is the start- 
ing point on the first circle C () ) plots the first segment. 

If there is a segment between two consecutive cir- 
cles in a path, this segment connects the ending point 
of the first circle with the starting point of the second 
circle. Therefore, automatically, the segments will be 
drawn when the circle path points are connected. 
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Lastly, connecting the final path point on the last 
circle C N _ , to the last point in a path plots the final 
segment P N . If there is a final segment P N beyond 
the last circle, then the coordinates of the last point on 
this segment are the coordinates of the center of the 
fictitious circle C N when R N — 0. Equation (13) 
gives the center of the fictitious circle C N as 


k N 


k N - 1 

n N - 1 

a N COS °N ~ 1 N sin °N 

k N 

> = < 

k N-\ 

> + r g < 

a N sin d N + l N cos 0 N 

. 0 . 


. 0 . 


0 


( 21 ) 


Temporarily, set R N = 0. This means that the tangent 
point of the segment P N with the fictitious circle C N 
is the same as the center of the fictitious circle. Also, 
this makes a N = R N _{, additionally, let e N = -1 so 
that 0 N = Ay . Then, the last point on P N is given by 


8. Applications in Path Design 

A computer program was written (in the C pro- 
gramming language) to compute some example paths. 
Two example runway exit paths are given for illustra- 
tion: (1) the first path is along two exterior tangent cir- 
cles and (2) the second approximates a spiral-type 
path. 

For simplicity, assume that the runway (think of 
P 0 in fig. 3) is parallel to the reference axis X ( . and 
points north (that is, T = y/ =0). Let the exit path 
be to the right (relative to motion just prior to the first 
circle). Although only examples of runway exit paths 
are given, the analysis applies equally well to a run- 
way entrance and takeoff path. 

8.1. Example 1: Runway Exit Path Using Two 
Exterior Tangent Circles 


(P 1 

v X last point 


h N - 1 
k N - 1 


> + R 


N - 1 
G 


0 


X 


C0S s ^ n ^7 V 

Rtf-l s * n + ^ N COS 


0 


( 22 ) 


Geometrically, equation (22) says that a vector from 
the reference frame to the last point in a path is the 
vector sum of a vector from the origin of the reference 
frame to the center of the last circle plus a vector 
(transformed to the reference frame) from the center of 
the fictitious last circle to the last point on the path. 
(See appendix C for specific equations to compute the 
beginning or ending point on any segment in a path.) 

7.3. Other Patterns Not Necessarily Smooth 


In reference 2, an airplane runway exit geometry 
has an initial constant radius of 2900 ft which makes a 
transition tangentially, after traveling 30°, into a final 
constant radius of 1500 ft. The travel angle on the sec- 
ond circle is also 30°. The final point on the exit path 
is approximately 2200 ft longitudinally and 600 ft lat- 
erally from the start of the exit turn. 

Example 1 is a path along two tangent circles 
(both exterior). Because the exit path is to the right, 
the direction of travel on the first circle is clockwise 
about Zq and the direction of travel on the second cir- 
cle is counterclockwise about Zj . The appropriate 
parameters to generate the path are as follows: 


Two circles, A 2 

Radius of first circle, Rq 2900 

Radius of second circle, R ^ 1500 

First circle is exterior, 1 

Second circle is also exterior, e ^ 1 

Angular travel on first circle, A^ , deg —30 

Angular travel on second circle, A-, , deg 30 


Smooth paths are the primary interest in this pub- 
lication, but the analysis applies to the design of a 
variety of patterns not just those that are smooth; see 
examples in appendix D. 


Assume that the exit turn starts at 3000 ft on the run- 
way. The resulting path, which agrees with that in ref- 
erence 2, is shown in figure 7. (For a left exit, set 
e Q = -1 and interchange the travel angles.) 
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Longitudinal distance along runway, x G , ft 


Figure 7. Example 1: two-circle runway exit path. 

8.2. Example 2: Approximating Spiral-Type 
Path 

Example 2 approximates the spiral-type path in 
reference 3. This approximation uses two circles at 
the beginning and at the end of the path with a linear 
segment in between. The first circle in each pair is 
exterior and the second is inward. The attraction of a 
spiral path is that it provides a smooth transition from 
the runway to the path. An airplane using a runway 
exit path that is initially spiral (decreasing radius) has 
a gradual buildup in the lateral acceleration for pas- 
senger comfort. Reference 3 uses 26 points (con- 
nected with linear segments) to approximate such a 
path to within a 2-ft error. Here, however, path param- 
eters were varied until there was a reasonable match 
between the approximate and the ideal paths. No 
attempt was made to design the path to within any par- 
ticular error criterion. The basic type of information is 
as follows (a user need not be aware of the particular 
variable symbol): 


Number of circles in path, N 4 

Exit direction, Exit direction = 1 Right 

Circle radius, ft, for — 

First circle, R Q 9100 

Second circle, R ^ 1700 

Third circle, R 1 1500 

Fourth circle, R^ 800 


Direction of motion on — 

First circle, (0 3 = -1 

Clockwise 

Second circle, to, = -1 

Clockwise 

Third circle, (0 3 = 1 

, Counterclockwise 

Fourth circle, (0 4 = 1 

, Counterclockwise 

Angular travel, deg, on — 

First circle, Aj = -4° 

4 

Second circle, A-, = -26° 

26 

Third circle, A 3 = 8° 

8 

Fourth circle, X 4 = 21 ° ... 

21 

Circle classification: 

First circle, e 0 = 1 

Exterior 

Second circle, e, = -1 .... 

Inward 

Third circle, c-, = 1 

Exterior 

Fourth circle, e 4 = -1 

Inward 


Segment length, ft, for — 

Initial segment (in front of first circle), / Q ..0 

Between first and second circle, / 1 0 

Between second and third circle, I 2 416 

Between third and fourth circle, / 3 0 

After last (fourth) circle, / 4 0 

The approximated path is shown in figure 8 and agrees 
well with the exact path. 

9. Following a Path 

A path is nothing more than a geometric curve 
(road). In trying to follow a path, one will deviate 
from the path; therefore, a correction is needed. A 
method to follow a path (refs. 1 and 2) is to orthogo- 
nally project the current position (moving point in ref- 
erence frame) onto the path as a target point and then 
make a correction to move toward this target point to 
reduce the lateral deviation from the path. (Think of 
tracking your shadow on the path which is always 
directly below or above you.) This approach appears 
successful in controlling lateral deviation and is the 
motivation for this section. However, actually control- 
ling motion to follow a path (refs. 1, 2, and 3) is not 
part of this publication. 

9.1. Target Point on Circle 

Imagine a new axis system (Z',T',Z') at the cen- 
ter of circle C, that is always parallel to the general 
reference frame (X ( .,Y (; ,Z ( .) as shown in figure 9. 
Let (-v G ,y G ,0) be a known, but arbitrary, point in the 
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800 |- 



general reference frame and project this point along 
the radius onto the circle. Refer to this point on the cir- 
cle as the target point. 





(-180° < f ; < 180°) (23) 


To compute the target point, first define a projection 
angle with respect to the (X'-,7'-,Z') axis system as 



Figure 9. Radial projection of known point to circle. 


The target point lies along the same radial line, from 
the center of the circle, as the known point, which 
means that the two points have the same projection 
angle. Therefore, in this axis system, the target-point 
coordinates are 


(x') t = R: sin (/)'■ 

' target ' r ' 

(24) 

(■y^) = R; COS <!)'■ 

''target ' r ' 

(25) 


and, in general reference frame coordinates, the target 
point coordinates are 

( x r •) = h: + (x') (26) 

v CM 'target ' v ''target v ' 

(y ) =*: + (/■) (27) 

v - Cr, /'target ' v - ''target v 
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where 


9.2. Target Point on Segment 

The starting and ending points on a segment 
define a line. This section gives a set of equations to 
compute a target point on this line. In general, if two 
straight lines intersect at right angles, the intersection 
point is the orthogonal projection of any point on 
either line onto the other line. Let ( x G j ,y G pO) and 
( x q ~>,y G 2>0) be two points that define a line L in the 
general reference frame (for example, the two points 
may be the beginning and ending points of a segment). 
Let be another line that (1) is perpendicular to L 
and (2) passes through an arbitrary but known point 
( X G fl’^G o’ 0)- The intersection of the two lines 
defines the orthogonal projection of the known point. 
The orthogonal projection is the target point. Obvi- 
ously, if L is parallel to X the target point is 


X G 

yc 

0 


target 




G,1 

2 G 

0 


( X G,1 “ X G, 2 ; 

y g = v G,a) 


(28) 


Otherwise, the target point is the intersection of the 
two lines: 


y G = mx G + (y GJ - my G l ) (29) 

2c = - l X G + (jo,a + \ X G,a J ( 30 > 

Call the intersection point the target point and make 
the notational change (x G ,y G ) = (x Ga ,y Ga ) to get 


■* G 

yo 

o 


target 


i 

. 2 

1 + m 


'"(yo-yo, t) + x o + m ~ x o, i 

X< , . 2 

m i x G + x o,\) + yo,\ + m yo 

o 


( x G,l* x G,l) (31) 


y G,2 ~ y G,1 

X G, 2 - X G, 1 


( 32 ) 


is the slope of the known line or segment and 
( x G ,y G ) now denotes the current position. 

9.3. Target Azimuth on Circle 

The azimuth of a moving target point on a circle is 
useful in correcting the heading of current motion. The 
azimuth at an arbitrary point on circle C, is 


Vi = V Xr , + V G i (i = 0, 1, AO (33) 

where y/ x is the known azimuth of the axis X G and 
y/ G i is the azimuth of the moving point relative to 
X G . Next consider how y/ (; ■ varies. 

The azimuth on circle Cj with respect to X (; is 
shown in figure 10(a) as a function of the projection 
angle </>'■ and the direction of motion on the circle. For 
example, if </>'• = 90°, figure 10(a) shows an azimuth 
of y/ G • = 90° for a clockwise motion, whereas 
figure 10(b) shows an azimuth of y/ (; • = 270° for a 
counterclockwise motion. 


Figure 10(c) shows a plot of \j/ G . as a function of 
(p'j and the direction of motion. Equations for com- 
puter computations are 


Vo,i = 180° ~ Q'i 
V G ,i = -ft 
Vo,i = 360° - <!>] 


(C. ccw; -180° < 0'. < 180°) 

(Cj cw; 0'<O°) ■ 
(Cj cw; <j>'j > 0°) 


(34) 

where b' is given by equation (23). The scenario is 
that current motion near a circle changes the circle 
projection angle 0' in equation (23), and this changes 
the azimuth in equations (34), which changes the tar- 
get azimuth in equation (33). 
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(a) Clockwise motion on circle. 


(b) Counterclockwise motion on circle. 



(c ) Azimuth angle as function of circle projection angle and direction of rotation on circle. 

Figure 10. Circle azimuth. 


9.4. Target Azimuth on Segment 

The constant target azimuth on a segment P • is 


W P . = ¥x g + tan 1 


%,2-VG,r 

x G 2 ~ X G, 1, 


(35) 


where ( x G j ,y G j ) is the starting point on the segment 
and ( x G 2 , >Y; t) is the ending point. 

Alternatively, the azimuth on a segment prior to 
the first circle is y/ p = y/ G + T ; the azimuth angle 
on any other segment following a circle is the same as 
the last azimuth angle on that circle (that is, the azi- 
muth at the ending point of the circle). 
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10. Switching Motion to Next Path 
Element 


next segment in the path P 1 and so forth. The loop 
continues until there are no more elements in the path. 


10.1. Sample Flowchart for Referencing 10.2. Switching Motion to Next Path Element 

Motion to Path 


Figure 1 1 shows a flowchart to reference motion 
to a path. The diagram sequences through the path 
along its various segments and circles. Assume there 
are N circles in a path, along with a final fictitious cir- 
cle (for convenience). Then, for each path index 
(/ = 0, 1, . .., N ) , there is a segment followed by a cir- 
cle. The fictitious circle follows the last segment. Ini- 
tially, the first element is P () , but if the length of this 
segment is zero, the first circle C Q then becomes the 
reference for the current motion. After completing 
motion on C () , the loop index increases by one to the 


Initial values: 



A user supplies some criteria to decide when to 
switch from one path element to the next. Numerous 
possible switching criteria (a research issue) may 
exist; for example, each element in the path has its 
own angle to go before the next element starts. When 
this angle reaches 0, motion switches to the next ele- 
ment and so on. 


This section defines an angle to go to the next suc- 
cessive element in a path. This switching angle is the 
central angle between the currently known position, a 
circle center, and the start (or end) of motion on the 
circle. Figure 12 indicates the switching angle for 
motion approaching an exterior circle. Bear in mind 
that the switching angle is only one indication of 
when to switch to the next path element; it does not 
indicate how close the motion is to the path at the 
instant of switching. Other parameters, such as the 
lateral position from the path, indicate the “closeness” 
of the motion to the path and also play an integral role 
in the switching decision. 

10.2.1. Switching Motion Away From Segment 

First, assume that motion is currently relative to 
the segment P Then, the angle to go tells when to 
switch motion to the upcoming circle C,. Let 
( x G ,y G ,0) be a known moving point in the general ref- 
erence frame. A transformation of this point to the cir- 
cle axis system is 


x i 


X G 

y,- 

,G 

■ = A r 

yc 

0 


0 

. i . 


i . 


(36) 


The radial projection angle of the moving point in the 
circle axis system is then 


0- = tan 




(-180° < </>• < 180°) (37) 


Figure 11. Sequencing motion along path elements. 


15 








Current location 



Figure 12. Angular location of current position from starting 
point on exterior circle. 

The central angle of travel before motion starts on cir- 
cle Cj is 

A,. = -(sign of <O /+1 0,.)( 1 +<? / )9O° + (y /+1 0 / (38) 

which is a function of the current projection angle, cir- 
cle classification, and the direction of motion. 
Equation (38) results by considering the various 
possibilities of approaching the start of motion on a 
circle. In essence, 

A • < 0 (Before motion starts on C-) 

A • = 0 (Start motion on C-) 

A - > 0 (Motion on C f ) 

(If there is a segment P N following the last circle in a 
path, a user can conveniently choose to set R N = j/^j 
to make A N vary from -45° to 0° as motion starts ana 
stops on this last segment.) Define an angle to go that 
is positive before motion is to start on C ■ as 

Angle to go = -A ■ (39) 

10.2.2. Switching Motion Away From Circle 

If motion is currently on circle C ( -, then the angle 
to go should indicate when to switch motion away 
from C • onto the next upcoming path element (circle 


or segment). In this case the angle to go is to the end 
of motion on C-. Although motion is on C -, A ( in 
equation (38) is positive and gives the central angle of 
travel on C •. Because the desired angular travel on C - 
is known, the angle to go for switching motion away 
from circle C ■ is 

Angle to go = |(A / ) max | - A,. (40) 

Normally, one knows the current position and then 
computes the angle to go to switch to the next element 
in a path. 

11. Comments on Control 

Knowing the angle to go and how this angle is 
currently changing with time, one can compute a time 
to go to switch to the next path element and, if desired, 
can start the switching maneuver a fraction of a sec- 
ond early for a better transition (less overshoot). In 
dealing with straight-line segments, reference 1 
divides the distance to go by the current speed to get a 
time to go to the end of the segment. When the time to 
go is within 0.05 sec, tracking switches to the next 
segment in the path. 

At a switching point in a path, there is a call for an 
instantaneous change in a velocity component; for 
example, in exiting from a runway onto a circular path 
element, there is a call for the velocity nonnal to the 
runway to change immediately from 0 to a value that 
is consistent with motion in a circle. Reference 2 
starts to build the nonnal velocity component slightly 
before reaching the exit for a smoother transition off 
the mnway, and reference 3 uses a spiral-type exit for 
smoothness. Error equations for controlling motion to 
follow a path are given in reference 1 . 

Reference 5 presents an algorithm, although not 
used, for lateral ride comfort limits. Current thinking 
(ref. 3) is that a lateral acceleration should not exceed 
0.1 5g and that lateral jerk (derivative of acceleration) 
should not be in excess of 0.05 g/ sec. 

12. Concluding Remarks 

This publication provides a general framework for 
the mathematical description of a smooth path along 
an arbitrary combination of circles and tangent lines. 


16 







The innovative approach is to model the path as a 
giant multilink manipulator and apply the mathemat- 
ics common to robotics. Any two successive circles in 
a path are either tangent or connected by a common 
tangent line. There may be an initial segment prior to 
the first circle in the path, and a path may begin or end 
on a circle or line. A user specifies a few parameters 
to configure the path. Necessary and sufficient 
conditions ensure that user inputs are consistent with a 


smooth path. An immediate application is designing 
smooth paths for use in automating the exit of an air- 
craft from a runway. Two example runway exit paths 
are given, and this publication defines an angle to go 
for switching from one path element to the next. 

NASA Langley Research Center 
Hampton, VA 23681-2199 
September 24, 1998 
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Appendix A 


Configuring Path With Circle Angles 

Appendix A illustrates geometrically how an initial configuration of circles and tangent lines changes with cir- 
cle angles in the design of a smooth path. However, in the main text circle angles are not explicit inputs in the 
design process. 

Figure A1 shows an example of configuring a path with circles and tangent lines. In figure Al(a), the initial 
configuration, prior to designing a path, is when the circle angles are zero. The configuration then changes with 
circle angles to design a path. The first two circles are exterior (the first by definition), and the third is an inward 
circle. There is a tangent segment prior to the first circle and another tangent segment between the second and third 
circles. Figure Al(b) shows the configuration after a rotation of = 1 10° about the center of the first circle, and 
figure Al(c) shows the final configuration after a rotation of 0 ~, = -110° about the center of the second circle. 
The example smooth path is then drawn along the circles and segments. The amount of travel along the last circle 
is 0 3 = 60°. 
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Linear 



(a) Zero circle angles. (b) Vary first circle angle. 



(c) Vary second and third circle angles. 

Figure Al. Constructing smooth path along circular arcs and tangent segments using circle angles. 
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Appendix B 


Homogeneous Transformation Matrices 

Homogeneous transformation matrices contain both orientation and displacement information between axis 
systems and are useful in following the path geometry. See reference 4 for application of homogeneous transfor- 
mation matrices to robotics. Consider the two axis systems a and b shown in sketch B1 : 



Sketch B1 

The two axis systems are related by the following 4x4 homogeneous transfonnation matrix: 


A b - 



1 

*5 -Si 

_ 

00 0 

1 


(Bl) 


where R b is a 3 X 3 rotational matrix that transforms a three-dimensional vector in axis system b to axis system a: 


and 


- 


r 


b 

X 

y 

r = RJ 

a 

y 

Z 


Z 


(B2) 


Pb = 


Px 

Py 

Pz 


(B3) 


is a 3 x 1 position vector directed from the origin of the axis system a to the origin of axis system b. 
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Bl. Inverse Transformation 


The homogeneous transformation matrix to transform quantities from axis system a to axis system b is the 
inverse of equation (Bl): 



(B4) 


where the superscript T represents the transpose of the matrix. 

B2. Transforming Points and Vectors 

A homogeneous transfonnation matrix transforms the coordinates of a point and the components of a vector as 
follows: 


' 


' 

X 



y 

■ =4-r 

v 

z 

z 

1 

1-1 

1 

s. 

^ J 


' 

V 


V 

V 


X 

V 

, i 

V 

y 

■ = A._ v 

? 

V 


V’ 

z 


z 

. 0 

i - 1 

0 , 


(Transforming coordinates of point) (B5) 


(Transfonning vector, such as velocity) (B6) 


B3. Specific Partitioning (Circle-to-Circle Matrix) 

The homogeneous transfomiation matrix from circle C ■ to the previous circle C •_ j is 


K-i j 

/-I 

Pi 

00 0 ! 

I 

1 


(B7) 


where 



cos 0 -sin 0. 0 

sin 9 cos 9j 0 
0 0 1 


(B8) 
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is an inverse Euler rotation matrix about Z ( , , and 


/-I 

Pi 


a i cos 0 . - /. sin 0 . 
a- sin 0 • + /• cos 0 • ' 
0 


is the position vector from the center of circle C _j 


to the center of circle 


C 



has 


(B9) 

components in 
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Appendix C 


Path Elements at Starting and Ending Points 

In switching from one path element to the next, one needs to know where the current element ends and the next 
begins. Appendix C gives equations for calculating the beginning and ending points on all circles and segments in 
a path. 

For a given path configuration, the homogeneous transformation matrix A' a from circle C to the reference 
frame is constant. In addition, for later convenience, define the matrix 



which is the same as matrix Aq 1 , except that / . +1 is 0 in A| +1 ; that is, except for the two elements: 


(Cl) 


^((U) = a i+l cos 0. +1 
l g + 1 ( I’ 3 ) = <* i+ 1 sin 0/+i 


(C2) 


Cl. Circles 

The starting point of motion on circle C ■ is 


x G,i 
* 

ycj 

0 

1 

where the superscript * identifies a starting point. If desired, equation (C3) expands to the more explicit vector expression: 



* 

x G,i 


Vi 

'+4' 1 ' 

cos Qj + Oj cos 0, -/,- sin 0 ? 

* 

yc,i 

= < 

k i- 1 

SjRj sin 0 . + a . sin 0 . + / . cos 0 ( . 

0 


0 

<. - 


0 


/q_l R i _ [ cos 0. - /. sin 0. 

k j _ l ' + R(; ' /? ( | sin 0 ( + l- cos 

oil 0 


(C4) 
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On the imaginary circle C N , the starting point is the same as the ending point because there is no motion on this 
circle; that is 


# 


* 

X G,N 


X G,N 

# 

» = . 

* 

yc,N 


- V G,W 

0 


0 


where the superscript # identifies the ending point. The ending points for all other circles are given by 


r 

# 


r v 

X G,i 


~ e i + l R i + l 

# 

y G,i 

T 

' _ l g 

o 

0 

i 


0 

1 

k 7 


(/ = 0, 1, ..., IV - 1) 


(C6) 


If desired, equation (C6) readily expands to 


' 

# 


x G,i 


# 

4 i 

y G,i 

’ = a g' 

0 


1 

< 



R i cos 1 

R i sin 9 i + 1 
0 
1 


(/' = 0, 1, IV - 1) 


(C7) 


C2. Segments 

The initial point on any initial segment prior to the first circle in a path is either known or computable from the 
exit coordinates, segment length, and azimuth angle. The final point on such a segment is the known exit point (or 
the starting point on the first circle). Hence, 

(/ , g) s t art = Known or computable point (C8) 


^o) end = Exit point (known) 


(C9) 


Thereafter, the starting point on segment P t is the ending point on the preceding circle: 




end 


# 

x G,i - 1 
# 

• V G,i— 1 
# 

z G,i-l 


(i'= 1,2, ...,N) 
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(CIO) 



The ending point on P- is the starting point on the upcoming circle: 


^)end = W 


start. 


x G,i 

% 

y G,i 


l Z G,i 


a = 1 , 2 , ...,n) 


(C11) 


As a matter of interest, after the first segment P Q , the relationship between the end point of segment P. 
(which is also the relationship between the ending point on C-_, and the starting point on C ) is 


(P) 

v ''start 


nd^G 1 


sin 0 ( . 
/■ cos Qj 
0 


(/= 1, 2, N) 


(Cl 2) 
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Appendix D 

Interesting Patterns, Not Necessarily Smooth 

One can generate some interesting patterns using circles and lines. Not necessarily designing patterns for 
smooth paths (as done in the main text) extends the variety of the patterns. The concern in appendix D is not 
whether a pattern is smooth but to generate some interesting patterns by varying mathematical parameters. 

For all patterns, except one, in this appendix, x exit = 3000, e 0 = 1 , and e N = -1 . Also, not shown, is the cir- 
cle direction of motion, which is the sign on maximum track angle; any other parameter not shown is 0. 
Equation (4) does not necessarily hold, which means there can be cusps in the pattern, but equation (3) still holds 
to give the signed segment length parameter. 


Dl. Patterns Using Four 90° Circular Arcs 

Each of the three patterns shown in figure Dl uses four 90° circular arcs. The first pattern (fig. Dl(a)) 
is a four-cusp pattern with N = 4, R 0 = R { = R 0 = = 2900, (^-i) max = (^ Li = - 90°, 

(X. 3 ) max = (A. 4 ) max = 90°, and e l = e 2 = e 3 = 1. The second pattern (fig. Dl(b)) shows a two-cusp pattern 




(c) Spiral. 

Figure D 1 . Patterns using 90° circular arcs. 


26 



resulting from shifting the lower portion of figure Dl(a) to the right by letting /-, = 6000. Figure Dl(c) is a spiral- 
ing pattern, with IV = 4, e, = e, = e, = -1, and (^.,) = (1 9 ) = (X-,) = (X A ) = -90°. 

° r 7 i z J x i 7 max v v 7 max v 


D2. Square- Wave Pattern 


Setting the radius of a circle in a path equal to zero directly connects its tangent segments. For example, with zero 
radii, the remaining parameters which generate the square wave are N = 5, e, = e 0 = e 3 = e, = 1, 


(X,) = (X 2 ) = (X 4 ) =-90°,(X.) =90°, /,=/,= 

v i^max v ^ y max v f^max _ x j'max 1 ^ 

the square -wave pattern, is shown in figure D2. 


1000, 


/ 3 = 2000 , 


and 


>4 = 


1000. This pattern, 


D3. Interconnecting Ring Pattern 

In the text, the circle classification parameter is either 1 or -1. Flowever, there is no need to impose this restric- 
tion in the quest for interesting patterns. For example, to draw two interconnecting rings, let = 0. (This 
removes all temis that multiply this parameter.) Other parameters for this pattern are N = 2, R 0 = R l = 1000, 
(/.j) = (A 2 ) = 360°, and 1 1 = 1000. The interconnecting ring pattern is shown in figure D3. 
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Figure D2. Square-wave pattern. 


Figure D3. Interconnecting ring pattern. 
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